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Abstract 


24th February 2015 


Focusing on the 2p — Is transition in atomic Hydrogen, we investigate through first order perturbation theory the time evolution 
of the survival probability of an electron initially taken to be in the excited (2p) state. We examine both the results yielded 
by the standard dipole approximation for the coupling between the atom and the electromagnetic field -for which we propose a 
c^ijf^ff-independent regularisation- and those yielded by the exact coupling function. In both cases, Fermi’s golden rule is shown 
tone an excellent approximation for the system at hand: we found its maximal deviation from the exact behaviour of the system 
t^e of order 10 ®/10 Our treatment also yields a rigorous prescription for the choice of the optimal cutoff frequency in the 
cljp^le approximation. With our cutoff, the predictions of the dipole approximation are almost indistinguishable at all times from 
t^^exact dynamics of the system. 

!^^^words: Spontaneous emission, Multipole coupling, Fermi’s golden rule, Zeno regime. 
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IN Introduction 

'description of spontaneous light emission by atomic elec- 
trg 2 ^ is a staple of textbooks in quantum physics [T] [2l |3] and 
qulmtum optics [4]. Oftentimes, several approximations are 
le in the treatment. The first one is the dipole approximation 
which consists [S] in considering that the decaying electron 
dcJSJ not emit light from its own position but rather from the 
p(5aJflon of the nucleus to which it is bound. This approxim¬ 
ation is usually justified by claiming that the electromagnetic 
wa^lengths which are relevant to the problem are much larger 
th^ the uncertainty on the electron’s position. It is thus argued 
t^jfe the electromagnetic field does not “see” the details of 
r tCTtj ter configuration at the atomic scale, and hence that the 
pt^ise location of the point of light emission is irrelevant. The 
s\s©nd approximation is often used to explain why the dipole 
c^^oximation can be made. It consists in noticing that, for large 
^qjigh times, the only field modes which effectively contribute 
t fTly )outaneous emission are the ones which are resonant with 
ttjp^ atomic transition frequency. This is known as Fermi’s 
goHlen rule. These resonant modes have a wavelength which is 
in^ed much larger than the relevant atomic dimensions, and 
th,eir interaction with the atom is very well described in the 
f^jij^ework of the dipole approximation. Hence to some extent 
tlfe-dipole approximation is justified by Fermi’s golden rule. But 
it^ a very general result [6] of quantum physics that Fermi’s 
linear decay cannot be valid at very short times. Indeed very 
short times obey what is called the Zeno dynamics, where the 
decay is always quadratic. One should therefore be skeptical 
of the validity of the dipole approximation at very short times, 
for which it may not be justified by the golden rule. Elsewhere 
[7] two of us have investigated the Hydrogen 2p — Is transition 
numerically, and developed a numerical method which enabled us 
to reproduce the dynamics of the system in the Zeno and Fermi 
regimes, but also at longer times in the Wigner-Weisskopf regime. 


In the following we investigate the validity of the dipole and Fermi 
approximations in the case of the 2p — Is transition in atomic Hy¬ 
drogen. In sect. we recall the main tools needed for the descrip¬ 
tion of spontaneous emission. In sect.[^ the central section of this 


manuscript, we show how Fermi’s golden rule emerges from rig¬ 
orous first order time-dependent perturbation theory, using two 
new independent arguments. First, in the dipole approximation. 



approximation and rigorously derive the short-time dynamics of 
the sytem. 


2 The decay of a two-level atom 
2.1 Position of the problem 

We consider a two-level atom, where the ground state | g) has 
angular frequency cug and the excited state state | e) has angular 
frequency uie, interacting with the electromagnetic field in the 
rotating wave approximation. We call e the (positive) elementary 
electric charge, me is the electron mass, and x is the position 
operator and p the linear momentum operator for the electron. 
The atom is considered to be in free space. The Hamiltonian H = 
Ha + Hr + Hi is a sum of three terms: the atom Hamiltonian Ha, 
the electromagnetic field Hamiltonian Hr, and the interaction 
Hamiltonian Hi. In the Schrodinger picture these read [8] 

Ha = Aug |g)(g| +huje |e)(e| 

(k)a(A) (k) 

A=±“' 

Hi = A (x, t = 0) • p 

me 

where A labels the polarisation of the electromagnetic field. In¬ 
troduce the polarisation vectors £(^=± 1 ) (k), which are any two 
(possibly complex) mutually orthogonal unit vectors taken in the 
plane orthogonal to the wave vector k. The vector potential 
A (x, t) is expanded (in the Coulomb gauge) over plane waves as 

A (x, t) = \j^ J [“(^) + h.c.] . 

( 2 ) 


(la) 

(lb) 

(l c) 
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Here 


dk= ^2n5{eo-k^)e{ko), (3) 

(27r) 

is the usual volume element on the light cone [9l [10] (where 6 
stands for the Heaviside distribution). Finally, the commutation 
relation between the photon ladder operators is given by 

[«(><) (k), aj;,) (q)] = 2||k||(27r)^5(k-q)(5^A. (4) 

The state of the system reads 

|V>(t))=Ce(t)e-‘“-‘ |e,0) 

+ /'<ifccg,A |g, lA.k) (5) 

A=±"' 

where | e, 0) means that the atom is in the excited state and the 
field contains no photons and | g, 1a, k) means that the atom is in 
the ground state and the field contains a photon of wave vector 
k and polarisation A. 


2.2 Importance of the interaction matrix element 


In the absence of any interaction (that is, if Hi were zero), the 
coefficients Ce and Cg,A (k, •) would not evolve and each term of 
the superposition Ht would oscillate at its own eigenfrequency. 
As such, the nontrivial features of the problem are encompassed 
by the matrix elements of the interaction Hamiltonian in the 
Hilbert (sub)space spanned by |e, 0) and |g, lA,k)- In the dipole 
^proximation, one approximates the imaginary exponential in 
M as e‘’‘'* ~ 1 (remember that it is the vector potential at 
t = 0 which appears in Hi, reflecting the fact that none of the 
Hamiltonians in Q depend explicitly on time). Then one is 
left to compute the relatively simple matrix elements of p. The 
result is well-known for the 2p — Is Hydrogen transition. Writing, 
for this transition, | g) =| Is) and | g) =| 2 pm 2 ), with m 2 the 
magnetic quantum number of the 2p sublevel considered, we have 
the following expressions for the electronic wave functions of the 
Is and 2 pm 2 sublevels: 



with ao the Bohr radius. The vectors ^^2 give the preferred dir¬ 
ectionality of the wave function of the 2p substates, the angular 
dependence of which is given by the usual spherical harmonics 
HD. They are given by 


^0 = Bz, 

(7a) 

^ Gx i IG-ii 

(7b) 


In the dipole approximation, the interaction matrix element then 
is [12] 


(Is, 1a, k |H/|2pm2,0) = — -(k) ■ im.2- ( 8 a) 

V eoc me do 3^ 

The dipole approximation is further discussed in sect. |3.1| For 
the same transition, if one keeps the full exponential in (|^, one 
obtains the exact matrix element m 


(Is, lA,k \Hi\2pm2,Q) = 


he 


22 


-(A) 


(k) ■ 


eoc meuo 3"^ 


[1 + (fao ||k||)^] 


(8b) 


Since we are interested in spontaneous emission, we set 
Ce (t = 0) — 1 and VA G {1, 2} Vk G K^Cg,A (k, t = 0) = 0. We 
want to compute quantities such as 

Pemiss.-»A,k (t) = |cg,A (k, t)|^ = | (g, lA,k | P (t) | e, 0) | (9a) 

where U (t) = exp |^(—i/fi) is the evolution operator for the 

system. It is well-known PE] that a time-dependent perturb¬ 
ative treatment of such a problem yields, to first order in time, 
+2 , 


■>-A,k 


W = I^S’sinc^ 


(c^o-c||k||) 
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(9b) 

where we introduced the notation ljq = tOg — cOg. From § we 
deduce that the probability that, at time t, the electron is still in 
the excited state, which we shall call the survival probability, is 
given by 


Psurv (t) = |Ce {t)\^ = ^ (k, t)| 


^ / dfe I (g, lA,k I Hr I e, 0) I 


A=± ■' 

(cuo -c||k||) ^ 

( 10 ) 

The integral in ( [10[ ) will be central for the rest of the article. 
Before we investigate it in further detail, let us remind the reader 
of the usual P 1131 [3] treatment of the problem. It consists in 
using the distributional limit 


lim ■ 

a —>- + oo 


(ax) 


= nS {x) 


( 11 ) 


of the square cardinal sine, to conclude that the only single¬ 
photon states available by spontaneous emission have a frequency 
equal to the atomic transition frequency ljJo- Indeed the use of 
in ([9b|l yields Fermi’s golden rule 


AA,k(t) ~ 27r-^ |(g, lA,k|H/|e,0)| 

t—>-+oo I I 


<5(a;o -c||k||). 

( 12 ) 

Notice that making use of (111 in (101 means taking the 
limit t —>■ -|-oo, a strange manipulation in the framework of 
time-dependent perturbation theory around t = 0. One then 
has to enter a subtle discussion of time regimes to conclude that 
Fermi’s golden rule is valid for times large enough to guarantee 
that one can approximate the square cardinal sine by its limit 
(111, but small enough to ensure that first-order perturbation 
theory still applies. This is discussed in more detail in [3] for 
instance. 


Our goal is not merely to point out the well-known fact that 
the golden rule is an approximation, but to acquire knowledge on 
what it is an approximation of. This is the question we investigate 
in the following sect.|^ where we focus on the case of the 2p— Is 
transition in atomic Hydrogen. 


3 Rigorous first order perturbation theory and 
Fermi’s golden rule 

We shall investigate corrections to Fermi’s golden rule from first 
order time-dependent perturbation theory. In sect. |3.l] we dis¬ 
cuss the validity of the ubiquitous dipole approximation for the 
atom-field coupling. In sect. [3.2] we propose a cutoff-independent 
regularisation procedure of the divergences which arise when the 
dipole approximation is made. We then go on to investigate the 
dynamics yielded by the exact coupling in sect. |3.3] 
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3.1 The dipole approximation: discussion 


It has been noticed that in the framework of the dipole approx¬ 
imation, it is necessary to introduce a cutoff over electromagnetic 
field frequencies. In our treatment, this can be seen as follows: 
write the survival probability of the electron in the excited state 
at time t as given by ( |8a[ | and (|10[ ). Performing the integration 
over the angles, we get, using (1^, the following expression: 

Psurv (t) = 1 - y doj a; sinc^ |^(cuo “ 2 J 


where a = e^/(47reo?ic) is the fine structure constant. The in¬ 
tegral on the right-hand side of (131 diverges at all times. In 
this situation, taking the limit t —>■ -l-oo under the integral to 
obtain Fermi’s golden rule is nonrigorous to say the least. Ac¬ 
cordingly, one introduces a cutoff frequency corresponding to an 
upper bound on the validity domain of the dipole approximation 
for the atom-field coupling. One then considers that electromag¬ 
netic held modes with frequency larger than this cutoff frequency 
are uncoupled to the atom. Several objections can be raised with 
regard to this procedure. We discuss them along with their more 
or less satisfying rebuttals: 




In order to introduce a cutoff, one must assume that high- 
frequency electromagnetic modes do not interact with the 
atom. This is justihed by a quick look at the exact coupling. 
The matrix elements are proportional to (Is | p exp (ik ■ x) | 
2pm2). When the wave number |jkj| becomes higher than 
the inverse of the Bohr radius ao, such matrix elements are 
de facto negligibly small because the oscillating exponential 
exp (ik ■ x) averages out during the integration, as seen from 


§. 


• The approximation yields cutoff-dependent results. Since 
the qualitative argument above only provides an order of 
magnitude (the ratio c/uq) for the cutoff frequency cue, this 
is especially problematic. Indeed, one can check that the 
truncated integral 


binary: the coupling function is taken to be exactly zero bey¬ 
ond the cutoff frequency. One could envision to introduce a 
smoother cutoff procedure along the lines of the cutting off 
of ultrarelativistic frequencies presented in [8], but we feel 
this would do nothing but introduce further arbitrariness in 
the model. 


3.2 The dipole approximation: regularisation 


For the reasons we gave in the previous sect. |3.1| we follow a 
different path: we refrain from introducing a cutoff and instead 
retain the result (131 from the dipole approximation without 
cutoff, but will endeavour to regularise the divergence in the 
integral. We will then extract the regular (finite) terms, and 
inspect them carefully. As far as we know this treatment of the 
present problem is novel. 


We focus our interest on the integral featured in (131, that is 

f-i-00 


/ 

J —< 


du) / (lj, t) where 


f = 6{uj)ujsmc \^uJo-uj)-j =6{uo)g{uj,t) . (15) 

Here 9 stands for the Heaviside step distribution. As we have 


discussed at length, (151 is a divergent integral. Nevertheless, 


we shall extract its regular (finite) terms. 


The idea here is that although the function does not 

belong to the vector space (R) of summable functions, it is a 
slowly growing function, and, as such, a tempered distribution. 
It therefore admits a Fourier transform in the sense of distribu¬ 
tions, which we write /(•,!). 


We thus compute the Fourier transform 




f 


da; / (cu, t) e 


(16) 


da; a; sinc^ ( (cjq — lj) ^ 


2 

¥ 


log ( —1 -I- — ) — [Ci ((a;c — a;o) t) — Ci (a;ot)] 

UJO 


LOo 


uJc — a;o 


(—1 -I- cos ((a;c — a;o) t)) + (—1 -|- cos (a;ot)) 


+t (Si ((a;c — a;o) t) + Si (oJot)) 


(14) 


and then shall take the limit r —>■ 0 at the end to retrieve 
the desired integral (15l. In this limit, some terms in / (r, t) 


become ill-defined, a consequence of the fact that is not 

summable. We will simply discard these terms at the end of our 
treatment, and focus on the well-defined terms in the limit r —>■ 0. 

The folding theorem and the well-known expression for the Four¬ 
ier transform of the Heaviside distribution yield, from (151 and 

(T^ 


where Ci stands for the cosine integral and Si for the sine 
integral [14] , is strongly dependent on the value of t he cutoff 
frequenc 3 jja;c. A discussion of the relevance of (14 1 is made 
in sect. HT In a recent paper [15] devoted to the decay of 
magnetic dipoles, similar questions were raised, and the 
Compton frequency was proposed as a cutoff. With this 
cutoff, Grimsmo et al. proposed a regularisation of their 
problem following the lines of Bethe’s mass renormalisation. 


/ {t, t) = 


5{-) - -vp- ) 

TT 


(t") 


(17) 


where vp stands for the Cauchy principal value of the subsequent 
func tion and the relation between g and g is the same as that 
(161 between / and /: 


9 (t t) = 




dcu g (w, t) e 


(18) 


From (15 I we get 


• When a cutoff is implemented, the distinction between elec¬ 
tromagnetic modes which are considered to be coupled to 
the atom and those who are excluded from the treatment is 

^We also note that when the interaction Hamiltonian is taken to be of 
the E • X form inst ead of the A • p form, the divergence of the integral cor¬ 
responding to i |l3[ l, which features instead of uj, is qua drat ic instead of 
logarithmic, and the dependence of the truncated integral 13 -which fea¬ 
tures the same substitution- on the cutoff frequency is very much enhanced. 
See [ 3 . 


9{T-,t) = 


/ +00 

dtJ- 

( 


(w — OJo) 


I g-iia{t+T)gi,aot _ 2e“ 


^). (19) 


We would like to compute this as a sum of three integrals corres¬ 
ponding to the three summands on the right-hand side of (191, 
but, taken individually, these integrals will diverge because of the 
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singularity at cu = ujo- The full integrand in (191, however, has no 


singularity at cu = oiq. Accordingly we introduce a small positive 
imaginary part e > 0 in the denominator, which enables us to 


compute (191 as a sum of three integrals. Introduce 


f 

J — c 


duj- 


{lo — iJo + ie) 


( 20 ) 


to rewrite 


g{T,t) = - ^ lim (e (r - t) + e‘"“*Ge (r + t) - 2Ge (r)) 

6 ^ 0 + V / 

( 21 ) 

Thus to compute g{-,t) we need only compute Ge, which we do 
now. Notice first that 

u) 1 oJo + ie 


(cu — oiQ + ie)^ a; — cuo + ie {oj — ojq + ie)^ 


Since e > 0, an application of Cauchy’s residue theorem therefore 
yields 

lim Ge (r) = 27re~'“°’^ (i + loot) 9 (r). (23) 


+0+ 


Plugging this back in (21 1 , we get 


g (r, t) = 27r^e [9 (r -t){i + uJo (r - t)) 

+9 (r + t) (i + LJo (r + t)) - 29 (r) (i + ojot)] . (24) 


Then we further plug this in (171 to get 
/’ + 00 


f{'r,t) = 


dcr - [6 (cr — t) (i + ojq (a — t)) 

T — (T 


(T,t)-^vp^ 

+9 (cr + t) (i + UJO (cr + t)) — 29 (cr) (i + tUQCr)] . (25) 

After some algebra we obtain 

= ^{ 7 r[ 6 '(T-t) (i + tJo (r-t)) 

+9 (t + t) (i + LOO (t + t)) 


—29 (t) (i + Wot)] — i vp 

/ /■* g-i^O^r rO g 

S-wo ( t J dcr-h J da- 




dcr- 


— da- 


f 

Jo 


-a— dcr 


(26) 


Hence we can finally write 


f {t = 0,t) = — {n \9 {-t) (i - Wot) + 9 {€) (i + wot) 


H-cjq I ^ 


— 2\6 ( 0 )] — i vp 

/: 


>(/: 


dcr- 


— da- 


dcr- 


- da 


+ / dcre 
0 


(27) 


Now we want to identify and discard the singular terms in \27\ . 
Making use of 9 {—t) + 0 (t) = 1, we can rewrite 20^) = 1 + 
9 (0) — 9 (—0). The quantity 9 (0) — 9 (—0) is singulaiQ and we 
discard it. We claim that the difference of twe two (principal 
value) integrals on the second line on the right-hand side of I 
also features a singular term. Let us show this. Rewrite 


/: 


dcr- 


L 


— da 


-<T Jo 'T 

-j-oo 


Now, for the n = 0 term in this sum, the integral diverges, and 
taking its principal value will not change that fact. Accordingly, 
we simply discard the n = 0 term in ( |28[ ) . We write the remainder 
of the series in closed form: 


+ 00 / -1 \n nt / 1 \n 


= Ci (wot) - log (wot) - 7 (29) 

where 7 is the Euler-Mascheroni constant. Computing the integ¬ 


rals on the third line of (271, we can rewrite the regular part of 
(27l as 


(22) /(T = 0,t) 


t 2 


-4sin^ ( ^ ) -I- 2 (Ci (wot) - log (wot) - 7 ) 


-l-TTWot ( sgn (t) -I-Si (wot) 


where r.p. stands for “regular part’. The term 


TT ^ ( sgn (t) -I- — Si (wot) 

t \ TT 


(30) 


(31) 


on the right-hand side of (301 is particularly interesting, and can 


be directly linked to Fermi’s golden rule. Indeed, remember from 
(131 that the d ecay probability (that is, 1 — Paurv (t)) features 
the product of (301 by t^. Further, notice that Si (wot) quickly 


converges to the Dirichlet value 7r/2 as t becomes substantially 
larger than 1/wo. For such times the leading term in (301 is 
clearly 27rwo/t (i.e., the limit of (311 as t —>■ -too) which is equal 
to the result obtained from illegally “sneaking” the limit t —>■ -|-oo 


into the divergent integral on the right-hand side of (131. Hence 


we have shown how the golden rule can be retrieved from a formal, 
cutoff-independent regularisation of the integral featured in the 
expression for the survival probability. We refrain from claiming 


that the terms on the first line of the right-hand side of (301 


are relevant descriptions of short-time deviations from the golden 
rule, as it is clear that using a more exact expression for the atom- 
field coupling will give better result^ This is examined in the 
upcoming sect. 3.3 Nevertheless, we can notice that (301 tends 
to zero as t —>■ 0 , which is an agreeable feature of our result. 

3.3 Exact coupling: vindication of Fermi’s golden 
rule 


In the present section we will start from the same integral (101 


and use the exact matrix element ( 8 bl. This treatment features 


no infinites and thus does not call for any regularisation proced¬ 
ure. Plus, it allows us to investigate short time deviations from 
Fermi’s golden rule in a more direct and reliable way. Start from 
the survival probability of the electron in the excited state 




^ + 00 

/ doj = 

Jo 




^ sine ( (wo - ^) 2 


(32) 


where we indroduced the notation wx = (3/2) (c/oq). The 
frequency wx is a natural cutoff frequency coming from the 


exact computation of the interaction matrix element ( 8 bL The 


o (-1)”, 2n 

n=0 ^ ’ 


f do-cr^"-!. 

integral in (32 1 

Jo 

^For the sake 


integral in (321 is finite at all times and we can compute it 


(28) 


^It can be understood as equal to sgn (0). 


result l |30[ l is plotted in Figs. and 
provide an accurate description of tn< 
system. 


where it is shown that it does not 
e very short-time behaviour of the 
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numerically or, as we shall now see, analytically. 


Cimoj 


Define 


/ + 00 

do; d (ui, t) where 

-OO 


d (cl), t) — 6 (uj) ^ ^ ^ sinc^ [ (coo — to) ^ 


It can be rewritten 

r + oo 


(33) 




do;- 



{uj - too) 




/ ^ —icjQt ^ — iojt^iiOQt 


- 2 ). (34) 


We wonld like to compute this as a sum of three integrals corres¬ 
ponding to the three summands on the right-hand side of (341, 


Figure 1: Jor dan loops in the complex Ai-plane used to compute the Fourier 
transform \3S\ . The poles — and iicjx of the integrand are represented 
by red circled crosses. 


where the at and bf coefficients depend on e. Whence the Fourier 
transform (|38|) 


but, taken individually, these integrals will diverge because of the 


singularity at cu = wo- The full integrand in (341, however, has no 


singularity at ca = aiQ. Accordingly we introduce a small positive 
imaginary part e > 0 in the denominator. Introduce 


H, it) = ^ 

/. 


ca 6 (ca) 

(ca — cao -b ie)^ 

1 + 

(it) 

2' 

4" 


he (t) = —2m 1^0 (t) e (ag + ai t) 

+d (t) 6-“^* (6j +btt + b^t^ + bt t^) 

-e {-t) e“>^‘ (bo -b br t -f b^ tV bj t®)] . (40) 


do; 


dtJ he (ca) 6 (ta) e 


(35a) 


(35b) 


One can then deduce He (t) and its limit as e —>■ O'*'. Setting 


+0+ 




bf 


+0+ 


Bf 


to rewrite 


(41a) 

(41b) 

(41c) 


Ip i^) ~ lini ^e (—t) -b e‘^°*He (t) — 2He (0)^ . one can see that 

(36) 

Thus to compute Ip we need only compute He, which we do 
now. The folding theorem and the well-known expression for the 
Fourier transform of the Heaviside distribution yield 


He it) = 


(5(-)- ;^vp- ) *hei-) 


it). 


We therefore need to compute the Fourier transform 

c+oo 


he it) = [ 

J —c 


dta he (ca) e 


(37) 


(38) 


B+ = Bg 

Bt = -Bf 

Bt = B^ 

Bt = -BT* = 


= Bo, 

= Bu 

= B2, 

= Bg. 


We give 


Ao = 


<^x (^x — 7tao) 


(cag+ca^) 


2 ^5 


Ai = -i- 


cao 


of he- We use Cauchy’s residue theorem. We know from ( |35[ ) 
that he has a second order pole at tug — ie and two fourth order 
poles at ±itax, pictured on Fig.[^ From ( |38[ ) we see that we have 
to close the integration path (Jordan loop) in the lower half of 
the complex plane for t > 0, and in the upper half of the plane 
for t < 0. 

It can be checked that the residues of he (ca) e““^* read 

Res (jie (•) e“‘'‘, cao — ie^ = (oq + ait), (39a) 

Res (^he (•) e“‘'f icax) = (bf -b b+ t -b bf -b bf t^) , 

(39b) 

Res (jie (■) e”’ *, —icax^ = (bg + bi t + b^ t^ + bj t^) 

(39c) 


1-b 


and 


Bg = - 
Bi = 


B2 = - 
Bs = 


(g;)_ 

cax (—Scag -b 30icaocax -b 48cax) 

96(cax -b icao)® 

cax (—3icaQ — 21caQcax -b 51icaocax -b 33cax) 
96(cax -b icao)® 

cax (—3icaQcax — 15caQcax -b 21icaocax -b 9cax) 
96(cax -b icao)® 

cax (—i^o^x — 3caQcax -b 3icaocax -b cax) 


(42a) 

(42b) 


96(cax -b icao)® 


(43a) 

(43b) 

(43c) 

(43d) 


The computation of Hg+ it) from (371 and (401 features no not¬ 
able conceptual or technical difficulty, but is quite tedious. We 
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only give the result, which reads 

Ho+ (t) = -ITT \d (t) e“‘““‘ (^0 + Ai t) 

+61 (t) e"“^‘ (B^ +Btt + Bt + Bt t^) 
-e i-t) e“>^‘ (Bo" +Bit + B^ + Bg- t^) ] 


+ e 


— iojQt 


—1 {Ao + Ai t) Si (cjot) H- A\ sin {ujot) 

(jjQ 


~ (^0 + Ai t) ("ci {iogt) + i—'(-j4ie'“°^ 

- [(BJ - B* t + Ba - B; t^) e““'^*Ei (uj^t) 

+ (Bo + Bi t + Ba + Bs e“^*Ei (-tJxt)] 

- — [(Bi + Bi*) - (Ba + Ba*) t + 2 (Bg + Bg*) t^] 

CJx 

-[(Ba — Ba) lo-xt + 2 (Bg — Bg) t]-(Bg + Bg) . 

CUx2 WxS 

(44) 


Here Ei stands for the exponential integral M- Now remains the 
unenviable task of adding three such terms as prescribed by (361, 


Ao + Bo + Bq — 0. 


(45) 


The link between (451 and the presence/absence of singular 


three poles), whence (45l. 


With the help of (45 I we can finally write, from (361 and (441 


If (i) = ^ |-2Ao ^log - Ci(a;o |t|)^ + itt (Bo - Bq) 


+Ai 


-4— sin^ ( ^ ) + iTvt ( sgn (t) + -Si (cuot) 
CJO V 2 / \ TT 


+i7r (t) [(^Bo* e'““* - Bo e"'““‘) 

+ (^-Bi + Bi t 

+ - Ba e"'“«*) + (-Bg* e‘““‘ + Bg e"'““*) t®) 

+e“^*6i (-t) ( (Bo - Bq 

- (-Bi* e“‘““‘ + Bi t 
+ (b *2 - Ba e‘““‘) - (-Bg* + Bg e'““‘) 

+e"“^‘Ei (a;xt) ((Bg* - Bg* t + Bg* - Bg* t^) 0^"“* 

+ (Bo - Bi t + Ba - Bg t®) 0 -'““^) 

+e‘"^*Ei (-a;xt) ((Bg* + Bg* t + Ba* t" + Bg* t®) e-“"“‘ 

+ (Bo + Bi t + Ba + Bg t®) 


( I Ce(f) 1-1 )M^ 



time(s) 


so as to obtain the exact expression for the survival probability 
as given by first-order time-dependent perturbation theory. Sim¬ 
plifications are frankly scarce here, but we are “saved” from the 
apparition of any singular terms by the identity 


ality I 

coupling as given by the perturbative solution j46[ (solid blue), and Facchi 
and Pascazio’s exact solution ITel (dashed green). The dimensionless con¬ 
stant — 6-4 X 10”^. Remember l/o^x — 1.18 X 10”^^ s. 


H—V [—2 ((Bi + B() oix — (Bg + Bg) dJx 

CJx 

+2(Bg + B3*))sin2 

+i ((Bg - Bg) iUx - (Bg - Bj)) uJxt sin (wot) 

+ (Bg + Bg) Wxi^ cos (wot) (46) 


terms is explained as follows. As seen from ( |39[ ), the quantity 
Ao + Bo + Bq is, up to a factor of ±2i7r, equal to the integral of 
ho+ over any closed curve V circling around the three poles of 
ho+ (see Fig. 0. We can take E to be a circle of radius R centred 
around z = 0. It is easy to see from Jordan’s lemmas that the 
integral of ho+ over such a curve vanishes when R —>■ +oc0 (and 
hence for any R large enough that the circle will still enclose the 


This expression, we are aware, is less than plain. Notice, how¬ 
ever, that it features the now familiar “seed” of Fermi’s golden 
rule, namely, the quantity 


If 2 

-iTrAi sgn(t)-|-Si (wot) 

t \ TT 


(47) 


found on the second line of (461. Notice from (421 that in the 
dipole limit wo/wx —>■ 0, one has Ai = —iwo, and one retrieves, 
from (321, the decay constant given by Fermi’s golden rule in 


the dipole approximation. Keeping wx to its actual value, we 
find a relative error |iAi — wo |/wq = 1.33 x 10“®. With wx 
and thus Ai kept to their actual value, our decay constant 
F = i 2^°/(3®7r) c^a^Ai matches that found by Facchi and 
Pascazio mnni, who treated the problem nonperturbatively, by 
finding the resolvent for the Hamiltonian 0 of the sytem in the 
Laplace domain and then getting back to the time domain. All 


the other terms in (461 are short-time deviations from Fermi’s 
golden rule. 


It is not too hard, but very tedious to show from (461, making 


use, of course, of (421 as well as (431, that the leading term in 


the Taylor series ofr If (t) at t = 0 takes the very simple form 
(wxt)^ /6. Whence the survival probability at very short times: 


(i) ~ 1 - 

t->o 3^7r 6 


(48) 


as deduced from (32 1 . This quadratic behaviour follows the lines 


of the usual Zeno regime [S] . The particular short-time expansion 
(481 has also been obtained by Facchi and Pascazio 11211161 . As 
shown on Figs. 0 and 0 the agreement between our perturbative 
treatment and the exact solution is very good. We can therefore 
use our method to investigate the short-time behaviour of this 
system in more detail. 


^For large uj, it is easy to see that hg+ behaves as lj 


The survival probability (32 I is plotted in Fig. 0 We see that 
the deviation from Fermi’s golden rule is of order 10“®/10“^. For 
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time(s) 


Figure 3: Behaviour of the survival probabi lity |ce(t)|^ with the exact 
coupling as given by the perturbative solution < |46^ {solid blue), and Facchi 
and Pascazio’s exact solution 1161 (dashed green). Here the exponential 
contribution Ze~^^ to the decay has been substracted, and only the non¬ 
exponential contribution thereto is plotted. The dimensionless constant 
^ (§)^ ^ 10“®. Remember l/cdx — 1-18 X 10“^® s. Also 

note Z ^ 1 - 4.39 A^ fT6l . 


1/CdX 1/9J0 



Figure 4: Decay of the survival probability |ce (t)|^ with the regularised 
dipole-ap pro ximated coupling (plotted for exhaustiveness) as given by ex- 

S ;sion l |30[ (dot-dashed green), the exact coupling as given by expression 
(solid blue) and Fermi’s golden rule (dashed black). The time axis is 
logarithmic. 


the sytem under study here, the golden rule is thus valid to an 
excellent approximation. 


To illustrate the transition between the Zeno regime (481 and 
Fermi’s golden rule as predicted by our treatment, we plot in 
Fig. the decay probability as given by ( |46[ ) as well as the short- 
time expansion (481 and the linear prediction of the golden rule. 
Note that the transition between the Zeno and Fermi regimes 
takes place around 10“^^ s after the start of the decay and that 
after s, the behaviour of the system is completely undistin- 

guishable from that predicted by the golden rule. 


4 Conclusion 

We have verified, in the framework of hrst order time-dependent 
perturbation theory, that Fermi’s rule that predicts a linear 
decay of the survival probability can indeed, for the 2p — Is 
transition in atomic Hydrogen, be called “golden”. The maximal 
deviation therefrom that we found is of order 10~®/10“^, a 
clear-cut endorsement. It is a success for a “rule” which, as 
was argued in sect. |2.2| is derived with the help of several, not 
obviously compatible conditions. As such, we think that much 
larger deviations from the golden rule could be found in other 


1/aJx 



Figure 5: Behaviour of the decay probability 1 — |ce (f)P with the regu¬ 
larised dipole- appr oximated coupling (plotted for exhaustiveness) as given 
by expre ssio n ( |30[ (dot-dashed green), the exact coupling as given by ex¬ 
pression l| 46t ( solid blue), Fermi’s golden rule (dashed black), and the Zeno 
behaviour Hsll (dotted red). Both axes are logatithmic. 


systems. As far as first-order perturbation theory goes, the 
crucial expression is the square cardinal sine integral in ( |10[ |. 
It tells us that with an enhanced coupling between an atom or 
another effective two-level system and (electromagnetic) modes 
which are off-resonant with the transition frequency of the 
two-level system, one would witness more important deviations 
from the golden rule. 


In [T^ Facchi and Pascazio also raised the question of the 
experimental observability of Zeno deviations from the golden 
rule. They considered the ratio between the Zeno time rz, where 
Pbutv (t) ~ 1 — (t/Tz)'^, and the lifetime l/F of the excited 

t—^■O 


level as the relevant parameter for the observability. We argue 
that the relevant ratio is that between the “cutoff time” rx 
and the Zeno time rz. The cutoff time is understood to be 
defined so that after rx, the system exits the Zeno regime in 
which the survival probability decays quadratically. Therefore, 
at t = rx, we have Paurv (t) = 1 — (rx/rz)^, and the strength of 
the Zeno effect is given by rx/rz. This is confirmed by looking 
at Fig. which shows that the maximal discrepancy between 
the predictions of the golden rule and the actual dynamics of the 
system is reached aproximately at the moment when the system 
exits the Zeno regime. There is a general method to obtain the 
Zeno time rz, which is centred [12] on the computation of the 
expectation value of the squared Hamiltonian of the system 
in the initial state. On the other hand, it is difficult to obtain 
the cutoff time rx without solving -at least perturbatively, as we 
did here- the dynamics of the system. Given the delicate nature 
of analytical approaches (see sect. |3.3| as well as (nmni) for this 
particular transition, which is much simpler than many other 
transitions to describe theoretically, we suspect that the best 
way to evaluate the cutoff time for a general transition would 
be by numerical evaluation of the integral on the right-hand 
side of (lOl. We showed in |7] that the ratio rx/rz scales, for 
hydrogen-like atoms with Z protons, like Z, a favourable scaling 
for the observability of the Zeno regime. 


Our investigations also shed light on the dipole approximation. 
We have seen that while the regularisation procedure of sect. |3.2| 
provides a nicely cutoff-independent treatment of the problem in 
the framework of the dipole approximation, and yields a result 
which is in agreement with Fermi’s golden rule at “long times”. 
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Figure 6: Behaviour of the decay pr obab ility 1 — |ce (t)P with the ex¬ 
act coupling as given by expression (solid blue) and the dipole- 

approximated coupling with cutoff frequency ujq — Wx/t/iT as given by 
expression (dot-dashed orange). Both axes are logatithmic. 


the predictions it yields on the very short time dynamics of the 
system are inadequate. Namely, it does not provide the correct 
dynamics in the Zeno regime, as seen on Fig. We might ask, 
however, how the predictions of the dipole approximation fare 
when the regularisation is performed more directly -and, argu¬ 
ably, less elegantly- via the introduction of a cutoff, as presented 
in sect. |3.l| As expected from the usual Zeno dynamics [B], the 
very short time behaviour yielded by the truncated integral (141 
is a quadratic decay in time. Namely 


1 ^ 2 3,2 


T 


dcj uj sinc^ 


(wo - w) 


2^° 2q3 2 


(49) 

This can either be seen by computing the Taylor series of the 
right-hand side of ( |14[ ), or, much more directly, by taking t — 0 
in the integral on the left-hand side of (49 b Now, remember 


that in the case of the exact coupling, the Zeno behaviour is 
given by (48 b One can then choose the cutoff frequency of the 


dipole approximation so that the very short time predictions of 
the dipole approximation, with cutoff, match the exact short time 


dynamics of the system. Comparison of (481 with (491 shows that 
a perfect match is reached if we choose 


wx c c 

wc = = —-~ .866—. 

V3 2 oo fflo 


(50) 


In Fig.j^we compare the predictions of the dipole approximation, 
with the carefully picked cutoff frequency (501, with the predic¬ 


tions obtained with the exact atom-field coupling. It is interest¬ 
ing, and quite impressive, that while the dipole approximation 
here was made to fit, by a simple choice of the cutoff frequency, 
the exact Zeno dynamics of the system, we see that with our 
choice for the cutoff, we obtain an excellent agreement during 
the transition between the Zeno and Fermi regime^ Hence we 
can conclude that for the 2p — Is transition in atomic Hydrogen, 
the dynamics of the system is very well described at all times 
within the framework of the dipole approximation, if one makes 


the “correct” choice (501 for the cutoff. 


^We also obtain an excellent agreement in the Fermi regime, but that 
was to be expected. The agreement is not perfect, though, as the decay 
constant in the e xac t and dipole coupling are slightly different. See the 
discussion below HD- 
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